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The measurements of CMB anisotropy have opened up a window for probing the global topology 
of the universe on length scales comparable to, and even beyond, the Hubble radius. For compact 
topologies, the two main effects on the CMB are: (1) the breaking of statistical isotropy in charac- 
teristic patterns determined by the photon geodesic structure of the manifold and (2) an infrared 
cutoff in the power spectrum of perturbations imposed by the finite spatial extent. We calculate the 
CMB anisotropy in compact hyperbolic universe models using the regularized method of images jlj|2| 
described in detail in paper-I |3j, including the line-of-sight "integrated Sachs- Wolfe" effect M, as 
Q\ • well as the last-scattering surface terms |1J. We calculate the Bayesian probabilities for a selection of 

ON \ models by confronting our theoretical pixel-pixel temperature correlation functions with the COBE- 

ON ■ dmr data. Our results demonstrate that strong constraints on compactness arise: if the universe is 

small compared to the 'horizon' size, correlations appear in the maps that are irreconcilable with the 
, observations. This conclusion is qualitatively insensitive to the matter content of the universe, in 

D ■ particular, the presence of a cosmological constant. If the universe is of comparable size to the 'hori- 

zon', the likelihood function is very dependent upon orientation of the manifold wrt the sky. While 
most orientations may be strongly ruled out, it sometimes happens that for a specific orientation the 
predicted correlation patterns are preferred over those for the conventional infinite models. The full 
Bayesian analysis we use is the most complete statistical test that can be done on the COBE maps, 
taking into account all possible signals and their variances in the theoretical skies, in particular the 
' high degree of anisotropic correlation that can exist. We also show that standard visual measures 

■^j- . for comparing theoretical predictions with the data such as the isotropized power spectrum Ce are 

not so useful in small compact spaces because of enhanced cosmic variance associated with the 
CNJ ' breakdown of statistical isotropy. 



I. INTRODUCTION 

The cosmic microwave background anisotropy is currently the most promising observational probe of the global 



spatial structure of the universe on length scales near to and even somewhat beyond the 'horizon' scale (~ cHq 1 ). As 
suggested by the concept of inflation, this relatively smooth Hubble volume that we observe is perhaps a tiny patch 
of an extremely inhomogeneous and complex spatial manifold. The complexity could involve non-trivial topology 
(multiple connectivity) on these ultra-large scales. Within a general program to address the observability of such a 
diverse global structure, a more well defined and tractable path would be to restrict oneself to spaces of uniform cur- 
vature (locally homogeneous and isotropic FRW models) but with non-trivial topology; in particular, compact spaces 
which have additional theoretical motivation [|] ||. For Euclidean (uniform zero curvature) or hyperbolic (uniform 
negative curvature) geometry, compactness necessarily implies non-trivial topology. Much recent astrophysical data 
suggest the cosmological density parameter in matter is subcritical ||, fl m < 1. Recently the presence of a significant 
cosmological constant (or, more generally, an exotic smooth component of matter) has been indicated by the high 
redshift supernova searches p0| and in studies combining large and intermediate angle CMB anistropy data with 



If this component does not compensate 
jolic spatial geometry for the universe; the 



observations of cluster abundances and large scale galaxy clustering [gd| 
for the deficit from unity, f2o = + < 1, this would imply a hyper: 
additional requirement of compactness then ushers into consideration topologically compact hyperbolic ( CH) universes, 
a field much richer in possibilities than the compact spaces with flat geometry when f2o is exactly unity. 

In a universe with non-trivial global spatial topology, the multiple connectivity of the space could lead to observable 
characteristic angular correlation patterns in the CMB anisotropy arising directly from multiple imaging of the source 
terms that give rise to the anisotropy in the CMB. Moreover, the modified structure of the eigenmodes in such 
spaces implies that angular correlations would differ from the predictions in the simply connected space with identical 
geometry (the latter plays the role of the universal cover of the multiply connected space), even in the absence of 
multiple imaging of the sources. In particular, compact universes cannot support modes whose characteristic length 
scale exceeds the linear size of the space; consequently, the inferred power of fluctuations in compact models at large 
scales would appear suppressed relative to the power on the universal covering space. A more subtle effect is the 
angular dependence of the theoretical temperature variance which reflects generic inhomogeneity in the topologically 
compact spaces. 
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Paper-I |3| describes the regularized method of images, a general technique that we developed [ljj2| for computing 
the spatial correlations in a universe with non-trivial topology. In this paper we shall apply the method to calculate 
the angular correlation of CMB anisotropics in CH universe models. The angular correlation between the CMB 
temperature fluctuations in two directions in the sky completely encodes the CMB anisotropy predictions of any 
model that postulates Gaussian primordial fluctuations. Using the full correlation function information on the COBE- 
dmr data, we have obtained limits on the size of flat torus universe models that are a factor of two sharper than 
obtained from the angular power spectrum alone ||. The main result was that the volume of the compact universe 
is constrained to be comparable to or larger than that of the observable universe. In jjj, we proposed on the basis of 
simple arguments that compact universes with hyperbolic geometry should be expected to respect similar constraints. 
We have now carried out the full Bayesian analysis on accurate CMB correlation predictions in a large selection of 
CH universes and can demonstrate that essential features of our general constraint is borne out. In this paper we 
describe the details of our CMB correlation computation, highlight the general correlation features in the compact 
spaces, describe the method of comparison to data and present the resulting constraints from COBE-dmr using the 
examples of a large and a small CH universe. A compilation of our constraints on the flat torus models and a large 



selection of CH models will be presented in a separate publication 13 1. 

The outline of this paper is as follows: In §[n] we recapitulate that the angular correlation function of CMB anisotropy 
on large angular scales, that predominantly arises through both surface and integrated Sachs- Wolfe effects, can be 
related to spatial correlations of the gravitational potential on the three-space hypersurface at the epoch of last 
scattering. This is a useful simplification in terms of computational costs for calculating the CMB correlation in 
compact spaces using our method. Although we restrict our calculations here to large angular scales, Appendix |A| 
discusses the implementation of the method of images to calculate the CMB anisotropies at smaller angular scales. 



In §111, we describe the computation of the angular correlation function of CMB anisotropy in CH models. The 
section also includes a quick review of some useful notions about compact spaces and the main result describing the 
regularized method of images from paper-I pj . 



In §IV, we discuss some of the typical correlation features in the CMB anisotropy that arise in compact universe 
models. We show that their origin is more readily understood by viewing the compact space as a tessellation of 
H 3 by the finite domains. In §fvl we present our results of full Bayesian probability analyses of large angle CMB 
anisotropy predictions for two CH models using the four year COBE-dmr data. In Appendix we show how to go 
beyond the Sachs- Wolfe effects to treat all aspects of CMB anisotropy using the method of images in compact spaces 
for high resolution CMB experiments. In Appendix [b|, we demonstrate that a byproduct of the statistical anisotropy 
of the CMB inherent in compact universe models is a considerably enhanced cosmic variance in the (isotropized) 
angular power spectrum Cg which completely characterizes the non-compact Gaussian models. This emphasizes that 
the pattern recognition aspect of the complete Bayesian testing of a model is essential to get the best constraints on 
allowed size of the compact space. 

II. CMB ANISOTROPY 

In the standard picture, the CMB that we observe is a Planckian distribution of relic photons which decoupled from 
matter at a redshift w 1100. These photons have freely propagated over a distance i? LS , comparable to the "horizon" 
size, a function of cosmological parameters. In a non-flat model the other length scale is curvature radius d c given 
by (c/ Hq)/ yjl — fio|. In a matter dominated (Qq = ^m) cosmology, i? LS ~ 2d c arctanh-^/I — f2o- For the adiabatic 
fluctuations we consider here, the dominant contribution to the anisotropy in the CMB temperature measured with 

1 /2 

wide-angle beams (6 FWHM comes from the cosmological metric perturbations through the Sachs- Wolfe effect. 

Although in this work we restrict our attention to large beam size and the Sachs- Wolfe effect, in Appendix |X] we show 
how effects which contribute to the CMB anisotropy at finer resolution can also be incorporated in the method of 
images. 

Adiabatic cosmological metric perturbations can be expressed in terms of a scalar gravitational potential <I>(x, r). 
The dynamical equation for ^(x, r) allows for separation of the spatial and temporal dependence in the linear regime,[] 
$(x, r) = (F(r) +£J(r))$(x), where $(x) is the field configuration on the three-hypersurface of constant time, whose 
amplitude is determined by the physics of the early universe. We use as our time variable dimensionless conformal 
time r, expressed in units of the curvature radius d c . Time dependence of the potential at the matter dominated stage 
is described by the 'growing' mode F(t) and the decaying mode E(t). In terms of the usual growth factor D{t) for 



1 At the scales appropriate to CMB anisotropies, damping effects on <3> can be neglected. 
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linear density perturbations, F = D/a, where a is the scale factor. The relative amplitude of the modes is determined 
by the matching condition at the moment t eq w 0.004/i~ 1 -v / 1 — ^o/^o of transition from extremely-relativistic to 
non-relativistic domination of the energy density. This gives E(t eq ) w F(t eq )/9. 

In this paper, we concentrate our study on open matter-dominated models with zero cosmological constant, fio = 
fi TO < 1, for which the growing mode evolves as E3] 



F(t) 



5 sinhr(sinhT — 3r) + 20(coshr — 1) 
(coshr — l) 3 



(1) 



in the matter dominated phase, r > r EQ . A non-zero cosmological constant can be trivially incorporated in our 
analysis by using the appropriate solution for F(r). 

We write the Sachs- Wolfe formula for the CMB temperature fluctuation, AT(q), in a direction q, in terms of the 
growing mode $ LS (x) of the potential at the three- hypersurface of constant time r = r LS , when the last scattering of 
CMB photons took place, $ LS (x) = F(r LS )$(x) 



dxf(x)$L S (qx) , fix) = 



F(r LS ) dr 



Fir) 



(2) 



where \ is the affine parameter along the photon path from X — at the observer position to x H = Rhs/d c - The first 
term is called the surface or "naive" Sachs- Wolfe effect (NSW). The second term, which is nonzero only if $ varies 
with time between t ls and now, is the integrated Sachs- Wolfe effect (ISW). The angular correlation between the CMB 
temperature fluctuations in two directions in the sky is then given by 

C{M') = (¥rti)^rtf)) = ^ LS (<zx H )<M9 x H )> 



\ T 



3./o 



dx fix) [(<M<7x)<I> L s(<7'Xh)) + (<M<7'x)<M<7X H )>] 



4/ d Xl fixJ I dx, fix,) (MgxJ^Ls^'xJ) 

lo Jo 



(3) 



The main point to be noted is that C(q,q') depends on the spatial two point correlation function, = 
($ LS (x)$ LS (x')) of $ on the three-hypersurface of last scattering. This is due to the fact that the equation of 
motion for <I> allows a separation of spatial and temporal dependence. As in eq.(||), the Sachs- Wolfe contribution to 
the correlation between temperature fluctuations in two directions in the sky can be split into three terms : 1) The 
surface term (NSW) which depends on the correlation between $ at the two points on the SLS; 2) the interference 
part correlating the value of $ at the points along one line of sight to the value at the SLS of the second line of 
sight; and, 3) an integral part which contains correlations between $ at points on one line of sight with those on the 
other. The last two terms constitute the ISW effect. If one considers the zero- lag correlation (the two lines of sight 
are identical), then the following holds: the first and third terms are positive definite, whereas the interference term 
comes in with a negative sign because dF(r)/dT is negative in the models that we consider here. 

In H 3 , the global isotropy of the space implies that the two point correlation function C(q,q') = C(9), where 
cos 9 — q ■ q' , and the CMB anisotropy can be described equally well in terms of its angular power spectrum d, 
defined by 



£ + 1/2 
£{£ + !) 



CtPt(cos6)W e 



(4) 



where the P# are Legendre polynomials and Wi encodes details of the experimental configuration, such as finite 
beamwidth. For COBE-dmr , Wt = B\ , where Bg is the beam, including a (sphericalized) approximation to finite 



2 A subtle aspect of the Sachs- Wolfe effect at large angular scales is that only the growing mode 7 ? (r L s)4 > (x) and not the total 
gravitational potential at last-scattering "^(x, r LS ) = (F(t ls ) + E(t ls ))&(h) contributes to the effect. The distinction is important 
for models with small value of £l m in which the time difference between the transition to matter domination at r EQ and the 
last scattering of photons at r LS is not large. With this important caveat, we shall still loosely call $ LS (x) the potential on the 
last-scattering hypersurface. 
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pixelization effects. In this paper, we use the experimentally-determined Wi for COBE-DMR , but, to set the scale, we 
note that a Gaussian Be ~ exp [—£(£+ l)cr^ cam /2] fit to the COBE-DMRbeam (including pixelization effects) gives 

The angular power spectrum of the CMB anisotropy on large angular scales arises mainly from the Sachs- Wolfe 
effect and is shown in Figure [l. The Ce contribution of the NSW term is plotted separately for comparison. The 
NSW contribution is suppressed at angular scales larger than the curvature scale due to focusing of geodesies in the 
hyperbolic geometry jllj and asymptotes to a constant value at small angular scales. The ISW contribution to Ce falls 
off roughly a factor of £ _1 faster than the NSW contribution. At small values of Slo, the ISW contribution dominates 
at large angular scales. 
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FIG. 1. The solid curves show the angular power spectrum Ce from the Sachs- Wolfe effect in the three infinite Ti. 3 universes 
with differing Qo. Also plotted as a dotted curve is the NSW contribution in each case. The curves are normalized such that the 
Ce from the NSW contribution goes to unity at large I. As fio increases towards unity, the relative ISW contribution diminishes 
and affects smaller values of I. While the ISW contribution is positive for Qq = 0.3 and Qo = 0.6, it is negative for fio = 0.8. 

Figure || illustrates some general features of the distribution of power in fc-space for the low multipoles of the Sachs- 
Wolfe CMB anisotropy, using the fourth multipole as an example. While the NSW contribution is always positive, 
the ISW contribution can be negative due to the interference term. The positive ISW contribution comes from larger 
values of k than for the NSW effect. Moreover, the small k contribution of the positive NSW effect is countered by 
the negative interference term in the ISW contribution (the second term in eq. (J3j) ) . Note the remarkable cancellation 
of the NSW contribution for f2o = 0.3. The relative contribution of the ISW effect decreases as f2o — ¥ 1- The presence 



3 When we speak about CMB anisotropy at "large angular scales" we always refer to the beam size (Tbeam and not the angular 
separation between the lines of sight, 8. This distinction is blurred in usual simply connected models. There, for example, 
the SW effect dominates both wide beam experiments such as COBE-DMRand also the correlation at large separation 9. In 
contrast, the distinction is important for compact spaces where points widely separated in angle may be physically close (see 
§^y|), resulting in the correlation between them being dominated not by SW but by short-distance effects, including the Doppler 
effect. It should also be noted that since compact universe models generically violate global isotropy, the decomposition (m) is 
not valid in those cases (the angular power in a multipole I is not evenly distributed in the azimuthal levels, m). 
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of the ISW contribution tends to relax the constraints on the size of a compact universe that can be directly inferred 
from the suppression of the low multipoles of the CMB anisotropy 16 1. 
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FIG. 2. The fc-space dependence of the integrand for the fourth multipole C4 from the Sachs-Wolfe effect in 3 infinite 
Ti? universes characterized by the values of fio shown. In each case the dotted and dashed curves are the NSW and ISW 
contributions, respectively. The ISW curve includes the NSW-ISW interference term as well. 



III. CMB ANISOTROPY IN COMPACT HYPERBOLIC MANIFOLDS 



A. Brief review of compact spaces 

We briefly recapitulate a few basic notions about compact universes that we discussed in paper-I ||. A compact 
cosmological model can be constructed by identifying points on the standard infinite flat or hyperbolic FRW spaces by 
the action of a suitable discrete subgroup of motions, F, of the full isometry group, G, of the FRW space. (The isometry 
group G is the group of motions which preserves the distances between points, i.e., leaves the metric unchanged). The 
infinite FRW spatial hypersurface is the universal cover, M u , tiled by copies of the compact space M. The compact 
space for a given location of the observer is most appropriately represented as the Dirichlet domain with the observer 
at its basepoint. Any point x of the compact space has an image Xj = 7^x in each copy of the Dirichlet domain on the 
universal cover, where ji £ T. The tiling of the universal cover with Dirichlet domains is a Voronoi tessellation, (a 
familiar concept in cosmology often used in modeling the large scale structure in the universe), with the seeds being 
the basepoint and its images. By construction a Dirichlet domain represents the compact space as a convex polyhedron 
with even number of faces identified pairwise under T. In cosmology, the Dirichlet domain constructed around the 
observer represents the universe as 'seen' by the observer and it proves useful in this context to define the outradius, 
i?> , the radii of the circumscribing sphere (smallest sphere around the observer which encloses the Dirichlet domain) 
and, the inradius, i?<, the radii of the inscribed sphere (largest sphere around the observer which can be enclosed 
within the Dirichlet domain) of the Dirichlet domain Q. Note that i?> and i?< are specific to the location of the 
observer within the compact space since the Dirichlet domains around different observers are not necessarily identical. 
An observer-independent (and Dirichlet-domain-independent) linear measure of the size of the compact space is given 
by the diameter of the space, dj^i = svLp x y£M d(x, y) fl7| , |l8|| , i.e., the maximum separation between two points in the 
compact space. 
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For cosmological CH models, A4 U = H 3 , the three-dimensional hyperbolic (uniform negative curvature) manifold. 
Ti 3 can be viewed as a hyperbolic section embedded in four dimensional flat Lorentzian space. The isometry group of 
Ti. 3 is the group of rotations in the four space - the proper Lorentz group, SO(3, 1). A CH manifold is then completely 
described by a discrete subgroup, T, of the proper Lorentz group, SO(3, 1). The Geometry Centre at the University 
of Minnesota has a large census of CH manifolds and public domain software SnapPea [|l9| . We have adapted this 
software to tile Ti. 3 under a given topology using a set of generators of T. The tiling routine uses the generator product 
method and ensures that all distinct tiles within a specified tiling radius are obtained. 

A CH manifold, A4, is characterized by a dimensionless number, Vm = YM/d 3 , where Vm is the volume of the 
space and d c is the curvature radius |20|. There are a countably infinite number of CH manifolds with no upper 
bound on Vm- The theoretical lower bound stands at Vm > 0.167 ]21j| . The smallest CH manifold discovered so far 
has Vm = 0-94 22 1. The Minnesota census lists several thousands of these manifolds with Vm u P to ~ 7. In the 
cosmological context, the physical size of the curvature radius d c is determined by the density parameter and the 
Hubble constant Ho: d c — (c/i?o)/Vl — ^o- The physical volume of the CH manifold with a given topology, i.e., a 
fixed value of VM/d 3 , is smaller for smaller values of Slo- 



B. Computing spatial correlation functions 

Here we summarize the main result of our regularized method of images that is discussed in detail in paper-I 
The correlation function on a compact space (and more generally, any non-simply connected space), M. = A4 U /T, 
can be expressed as a regularized sum over the correlation function on its universal cover, M. u , calculated between x 
and the images -fx' (7 G T) of x': 



a(x,x')=^a(*,7x') 



7er 



= ^«(x,7x')--L / dx'Q(x,x' 
Vm Jm- 



(5) 



The local isotropy and homogeneity of Ti. 3 implies £$( x ? x') depends only on the proper distance, r = d(x, x'), 
between the points x and x' . The eigenfunctions on the universal cover are of course well known for all homogeneous 
and isotropic models 23 1. Consequently £$(x, x') can be obtained using its eigenmode expansion. The initial power 
spectrum P$(fc) is believed to be dictated by an early universe scenario for the generation of primordial perturbations. 
We assume that the initial perturbations are generated by quantum vacuum fluctuations during inflation. This leads 
to 



a(x,x') = sw = 



df3 (3 sin(/3r) 



(/3 2 + 1) /3sinhr 



(6) 



where (3 = y/(kd c ) 2 ~ 1 and P$(/3) = (3{f3 2 + l)P<s,(k)/(2ir 2 ). 

In the simplest inflationary models, V^, the power per logarithmic interval of k is approximately constant in the 
subcurvature sector, defined by kd c > 1. This is the generalization of the Harrison-Zeldovich spectrum in spatially flat 
models to hyperbolic spaces p^ , Pq ]. Sub-horizon vacuum fluctuations during inflation are not expected to generate 
supercurvature modes, those with kd c < 1, which is why they are not included in eq.(^). Indeed, since H 2 > l/(ad c ) 2 , 
for modes with kd c < 1 we always have k/{aH) < 1 so inflation by itself does not provide a causal mechanism for 
their excitation. Moreover, the lowest non-zero eigenvalue in compact spaces, k\ > 0, provides an infra-red cutoff in 
the spectrum which can be large enough in many CH spaces to exclude the supercurvature sector entirely (k\d c > 1). 
(See [|.) Even if the space does support supercurvature modes, some physical mechanism needs to be invoked to 
excite them, e.g. ; as a byproduct of the creation of the compact space itself, but which could be accompanied by 
complex nonperturbative structure as well. To have quantitative predictions for P$(fc) would require addressing this 
possibility in a full quantum cosmological context. We note that our main conclusions regarding peculiar correlation 
features in the CMB anisotropy (see § IV) would qualitatively hold even in the presence of supercurvature modes. 

Although equation (^) encodes the basic formula for calculating the correlation function, it is not numerically 
implementable as is. Both the sum and the integral in equation (^) are divergent and the difference needs to be taken 
as a limiting process of summation of images and integration up to a finite distance r* : 



= lim 



E 



Air 



dr sinh 2 r ^$(r) 



d(x, 7i x') 



(7) 



G 



The volume element in the integral is the one appropriate for H 3 . Numerically we have found it suffices to evaluate 
the above expression up to r* about 4 to 5 times the domain size i?> to obtain a convergent result for £J(x, x'). 



C. Computing the CMB correlation 

For Gaussian perturbations, the angular correlation function, C(q,q'), completely encodes the CMB anisotropy 
predictions of a model. To make maps, the celestial sphere is discretized into N p pixels labeled by p. N p is determined 
by the angular resolution, ~ 7° degrees for COBE-DMR. The COBE-DMRmaps had (2.6°) 2 pixels, corresponding to 
N p = 6144, though compressing the data into N p = 1536 pixels leads to no information loss. We use this number 
of pixels in § [y], where we confront our models with the COBE-DMR observations using Bayesian statistical methods, 
and also in the maps shown in this paper. 

For Gaussian statistics, the pixelized theoretical maps are fully determined by the N p x N p pixel-pixel correlation 
matrix, Ct pp > = C(q p , q p >). The expression for C(q, q') in eq. (|^) involves an integral along the line of sight from the 
observer to the surface of last scattering - the integrated Sachs- Wolfe term (ISW). We find that a simple integration 
rule using Nl points along each line of sight gives accurate results when the points are spaced at equal increments AF 
of F(t). Consequently, to get the full Ct PP ' matrix we need to evaluate the correlation function, £|(x, x'), between 
N p (N p + l)Nf j /2 pairs of points on the constant-time hypersurface of last scattering: 

1 2 AF Nl fAF\ 2 Nl Nl 

C TpP ' = ^UXhQp, Xn%') + 3 — K*(Xh9p,X»9p') + tUXiQp' , XhQp)] + 4 ( — J ^ ^ Wij^UXitip', Xj<h>) 

i—1 ^ ' i—1 j — 1 

(8) 

where Wi and Wij are 0(1) coefficients that depend on the specific difference scheme used. 

To calculate Ct pp ' in the standard infinite open cosmological models with this real-space integration to an accuracy 
comparable to that of the traditional evaluation in fc-space, Nl ~ 10 is sufficient. Moreover, the real space integration 
is much faster in terms of CPU time. For the method to remain accurate in compact models, Nl should, of course, 
exceed the number of the times a typical photon path crosses the compact Universe. We found Nl ~ 10 is still enough 
for the models we have analyzed so far. 



IV. CORRELATION FEATURES IN THE CMB ANISOTROPY 



In standard cosmological models based on a topologically trivial space such as Tt 3 , the observed CMB photons have 
propagated along radial geodesies from a 2-sphere of radius R LS (that we refer to as the sphere of last scattering, SLS) 
centered on the observer. The same picture also applies to CH models when the space is viewed as a tessellation 
of the universal cover tiled by the Dirichlet domain with the observer at the basepoint. If Vls, the volume of the 
SLS, is much larger that that of the compact space, the photons propagate through a lattice of identical domains. 
As a consequence, strong correlations build up between CMB temperature fluctuations observed in widely separated 
directions. The correlation function C(q,q') is anisotropic and contains characteristic patterns determined by the 
photon geodesic structure of the compact manifold. These correlations persist even in CH models whose Dirichlet 
domains are comparable to or slightly bigger than the SLS. This is the key difference from the standard models where 
C(q, q 1 ) depends only on the angle between q and q' and generally falls off with angular separation. 

The pattern of strong correlations is directly linked to the way the points (on the universal cover A4 U ) enclosed 
by the SLS are equivalent under the topological identification T. Consider a set of points S on the universal cover. 
Then for all 7 £ T the pair of subsets 7^ n S and S n 7~ 1 S' are equivalent under T, i.e., either contain the same 
points of the compact space M. or they are both empty. (Obviously, the isometry 7~ 1 G T maps the first set onto 
the second.) If S consists of the points on the SLS, a non-empty set S D 7^ would be a circle on the SLS which is 
pointwise identical to another circle S f) J~ 1 S on the SLS. If the CMB anisotropy is dominated by the surface term, 
then the CMB temperature along one circle is expected to be identical to that along the other circle. This effect 
was first understood and much emphasized in the works byCornish, Starkman and Spergel |^8| . The necessary and 
sufficient condition for the existence of matched circles (non empty subsets S (1 jS) on the SLS is that R hs > i?< , 
i.e., the SLS is not completely enclosed within the Dirichlet domain. We discuss the CMB correlations in our models 



that arise from this kind of identification in S IV A . 
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FIG. 3. The origin of CMB correlation patterns in a multiply connected universe is illustrated. The compact space is 
represented by the (hexagonal) Dirichlet domain, D, around the observer, O, which tiles the universal cover. The tiles 7D and 
7 _1 D under one of the face translations 7 6 F and its inverse are shown. The points O' and O" are the images of the observer, 
O, under 7 and 7 -1 , respectively. The point B' is the 7-translate of B and the point A' is the 7 _1 translate of B. The solid 
circle Sls represents the sphere of last scattering (SLS) on the universal cover; the dashed circles, labeled jSls and 7 _1 Sls i 
are the translates of SLS under 7 and 7 -1 , respectively. The intersections of the sphere Sls with ^Sls and ~/~ 1 Sls create a 
matched pair of circles on the SLS. Similarly, the shaded lens-shaped regions within the SLS consist of identical sets of points. 
Consider the lines of sight OA and OB directed towards the centers of the matched circles. The ray O'B' is the image of OB 
under 7 and the ray O" A' is the image of OA under 7 -1 . By the arguments outlined in the text, the segment AB' = OAnOB' 
is identical to A'B = OB n OA'. 



The CMB anisotropy has contributions from the entire line of sight and it is useful to study the subsets of identical 
points enclosed by the SLS. If S is the set of all points enclosed by the SLS, then the pairs of subsets of S that are 
identical are lens-shaped regions created by the intersection of two balls (See Figure ||). Again the necessary and 
sufficient condition for the existence of such identified regions within the SLS is R LS > R < . Further if the SLS does 
not enclose more than one Dirichlet domain, these lens-like regions are in the direction of the faces and are adjacent 
to the SLS. (We find from our analysis of the COBE-DMR constraints that viable CH models have volumes comparable 
to or more than that enclosed within the SLS. See § |v|.) In these cases it is clear that correlation patterns are built 
from points close to the SLS and directly reflect the shape of polyhedral Dirichlet domain. 

The ISW-CMB correlation depends on the correlations between the points lying on the two (radial) lines of sight 
from the observer to the SLS. Hence it is instructive to study whether two different lines of sight (radial lines) in the 
SLS contain a set of identical points. Consider the set of points L\ and L2 along two distinct lines of sight. Then for 
all 7 € F, the subset L\ [I7 £2 of L\, if non-empty, would contain the same points as the subset 7L1 n £2 of L^. The 
necessary and sufficient condition for the existence of such pairs of lines of sight is again R hs > i?< . It is possible to 
show that there are infinite pairs of lines of sights which share at least one common point, but what is more interesting 
and relevant are the pairs where the identical subsets contain a segment. It is a straightforward exercise to verify 
that every pair of lines of sight pointing towards centers of matched circle s discussed above must contain a segment 
consisting of identical points (See Figure ||). As we shall discuss in § [VB, this result is important in understanding 
the correlation patterns in CMB anisotropy when the integrated Sachs Wolfe contribution is significant. 

In the regime where R LS < i?<, there are no points within the SLS that are topologically equivalent and the 
specific correlation features discussed above are absent. However, quantitatively the CMB correlations continue to 
have observable deviations from that expected in a simply connected model until R LS is substantially smaller than 
R < . Typically the correlation pattern around any point in the sky is anisotropic and distorted. 

In simply connected models of the universe, the global homogeneity and isotropy ensures that the CMB sky is 
statistically equivalent for all observers. These global symmetries are generically absent in compact (open or flat) 
universe models. Consider two distinct observers. If one of them measures the correlation function between pairs 
(x, x'), the corresponding measurement for the other one will be between pairs (gx, gx.') where g is the element of the 
isometry group of the universal cover, g £ G u , which transports the first observer to the position of the second one. 
(The motion g always exists since the universal cover is homogeneous.) For example if x, x' belong to the SLS as seen 
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by first observer, then gx, gx' are on the SLS from the second observer's point of view. Each pairwise correlation 
value £|(x, x') ^ s determined by the set of distances T>(x, x') = {d(x, 7jx'), 7j £ T}. Moving the observer translates 
this set to {d(<7X, 7j-gx/)}. Since G u is the isometry group on the universal cover, it conserves distances; in particular, 
d(gx, 7jgx/) = d(g~ 1 gx, g~ 1 jjgxf). Thus, under a general relocation of the observer by g, the set of distances £>(x, x/) 
is transformed as {rf(x, 7jX/)} — > {d(x, g^^trx/)}. 

If r is a normal/invariant subgroup of G u , (g _1 7<? £ L,V7 € T, g G G u ), moving the observer simply reshuffles the 
terms in the set T>, and £|,(x, x ') = £!>(5 x j 5 x ')- This is the case with simple tori, in which opposite sides are identified 
by pure translation without rotation. In the case of CH spaces (as well as tori with twists), T is not a normal subgroup 
of G u . Equivalent only are those observers mapped by a g belonging to the isometry group G of the compact space 
M. itself, not in the larger isometry group of its cover M. u . Each element of the isometry group G on M. commutes 
with all the elements of T, with the immediate result that the Ct pp ' calculated will be the same up to rotations of 
the sky if the observer is moved along the orbits of G in M.. 

In§" 
space. 



IV C, we discuss the correlation pattern in the CMB that arises due to the global inhomogeneity of the compact 
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FIG. 4. This shows the behavior of the correlation function of the CMB temperature on a great circle in the sky in a CH 
model (m004(-5,l)). The solid curves in the first and third panels show C{6) for the surface term in the Sachs-Wolfe effect 
(NSW) for fio = 0.3 and 0.6, respectively. They reflect the spatial correlation along a circle on the sphere of last scattering 
(SLS). The number of peaks in C(8) matches the number of Dirichlet domains that the circle intersects. The smooth curves 
show corresponding results for the simply connected infinite TC 3 models with flo = 0.3 and 0.6. The second and fourth panels are 
analogous to the first and third panels, except with the integrated Sachs- Wolfe (ISW) effect included. In this more physically 
correct case, the sharp NSW peaks have been diluted by the ISW contributions. However, the ISW effect induces new features; 
in particular, note the appearance of strong negative correlations in panels two and four. 



A. Correlations due to the NSW Surface term alone 

In Figure U the complex behavior of the correlation function is illustrated with the example of the small CH model 
(m004(-5,l)). The SLS encompasses 150 domains for J7 = 0.3 and « 20 domains for O = 0.6. It is comparable 
to the size of one domain for Qq — 0.9. The angular correlation along any arbitrary great circle in the sky in lower 
f2o models shows distinct peaks as one encounters repeated copies of the Dirichlet domain. The peaks are more 
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pronounced when one considers only the surface terms of eq.(|2|) in the CMB anisotropy. Including the line of sig ht 
ISW contribution tends to smearing the peaks, but also adds its own characteristic features as discussed in § IV B . 
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FIG. 5. This shows two full-sky (noiseless) CMB anisotropy maps, plotted as pairs of 180° diameter hemispherical caps, 
one centered on the South Galactic Pole (SGP) and one on the North (NGP). They are one of an infinite number of possible 
random realizations based on the computed pixel-pixel correlation matrix for the model in question. Both surface and integrated 
Sachs- Wolfe effects have been included. The power was normalized to best match the COBE-DMRdata. The contours are linearly 
spaced in 30 fiK steps. In contrast to Figs. [u] and [ll], the maps are not optimally filtered. The model labels L(arge)CH and 
S(mall)CH refer to the CH models v3543(2,3) and m004(-5,l), respectively. (The model number associated with the topology 
corresponds to that of the census of CH spaces from the Geometry center, Univ. of Minnesota; SCH is one of the smallest and 
LCH is one of the largest spaces in the census.) The value of Qo in each was chosen so that Vm ~ V^ LS . The matched pairs of 
circles expected if the CMB anisotropy is dominated by the surface terms [28] are superimposed on the map for each model. 
Each pair is labeled by the same number centered on the circles. The relative phase is shown by identified points marked by a 
diamond and a triangle on each circle in a pair. For clarity, we show only the eight largest pairs out of 35 for the LCHn 0= o.6 
case. Even at COBE-dmr resolution, we find the cross-correlation between the temperature along matched circles is very good 
in the SGHn 0= o.g model. The ISW contribution is larger at Qo = 0.6, and the cross-correlation coefficients are systematically 
smaller for LCHn =o.6 circle pairs. 

Consider compact universe models which are small enough such that the SLS does not completely fit inside one 
domain, R hs > R < . The CMB temperature is expected to be identical along pairs of circles if temperature fluctuations 
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are dominated by the surface terms at the SLS |28j . We identify these matched circles on the sky in our models and 
check the extent of cross correlation seen at the angular resolution of COBE-dmr . Figure [| shows the matched pairs of 
circles in two ( "small" and "large" ) CH models superimposed on random realizations of the theoretical sky generated 
from the appropriate pixel-pixel correlation matrices. The models were chosen to have a volume comparable to the 
volume Vls within the SLS. Even with the coarse pixelization of COBE-DMR, we do see fairly good cross-correlation 
in the CMB temperature along matched circles in our realizations. Again, the pattern of correlated circles is more 
pronounced when the Sachs- Wolfe surface term is dominant, as happens in the small model (SCH) with f2o — 0.9. 
We compute the cross-correlation coefficient between the temperature fluctuations along two matched circles C\ and 
C 2 = iC u 

Pl2 = (AT(x)AT( 7 x))/ [<AT(x) 2 )(AT( 7 x) 2 )] 1/2 , x e d, 7 x e C 2 . (9) 

(In the case of a single realization we replace the statistical average by the integration over x along the circles.) In 
the SCH £lo = 0.9 case, is in the range 0.6 — 0.95, whereas it is in the range 0.2 — 0.6 for the large model (LCH) 
with Qq = 0.6. 

The circles of identified pixels on the SLS are not the whole story. Enhanced NSW cross-correlations, but at a lower 
level, also exist between all pairs of points on the SLS which are projected close to each other on the CH manifold. This 
can be seen as secondary maxima in the examples in Figure^ these are absent in standard cosmological models. Some 
features persist at a detectable level even when the compact universe encompasses the SLS, i.e., R LS < R < , although 
circles are absent in this case (the effect dies out for R hs <C R<). When the relative ISW component is significant, the 
geometrical patterns based on pointwise identifications on the SLS are supplanted by more complex features arising 
from identifications between photon geodesies, e.g., the strong negative correlations evident in FigureS. 
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FIG. 6. In the eight panels of the figure, the CMB temperature along the 8 matched pairs of circles of the SCH model 
shown in Figure ^ (top panel) is plotted as a function of a polar angle. The temperature shown is from the same realization 
of the CMB sky that is used in Figure [| and the diamond and triangles are in correspondence. The CMB anisotropy includes 
both the surface and integrated Sachs- Wolfe effects. At fio = 0.9, the ISW contribution is small, and thus the values of the 
correlation coefficient are high. 
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FIG. 7. In the eight panels of the figure, the CMB temperature along the largest 8 matched pairs of circles of the LCH 
model (fio = 0.6) shown in Figure ^| (bottom panel) is plotted as a function of a polar angle. The temperature shown is from 
the same realization of the CMB sky that is used in Figure |H| and the diamond and triangles markers are in correspondence. The 
CMB anisotropy includes both the surface and integrated Sachs- Wolfe effects. At f2o = 0.6, the ISW contribution is significant, 
resulting in small values of the correlation coefficient. 



B. Correlations with the ISW effect included 



In a hyperbolic model, the temperature fluctuations depends on the entire path of the photons from the sphere of 
last scattering to the observer. The relative contribution of the ISW effect increases as the value Sl decreases. As the 
ISW component increases, it tends to wash out the patterns arising solely from the NSW (described in the previous 
section). However, it also introduces additional features in the CMB correlation. In Figures || and [?], the temperature 
fluctuations along matched pairs of circles in a realization are plotted. Since Qq is closer to unity in the small model 
(SCH), the temperature fluctuations along the circles are more tightly correlated than in the large model (LCH). 

An example of ISW-induced features is the significant negative correlations between widely separated directions 
in the sky, as is clearly seen in the panels two and four of Figur e |} In fact, the highly anti-correlated regions tend 
to lie at the centers of the matched circles discussed in § IV A[ Figure || plots the anticorrelated pairs of pixels 
separated by more that ss 10° below a threshold in the correlation coefficient of C vv , / ^/C vv C p , v i = —0.3 in the LCH 
model with f2o = 0.6. This should be contrasted with the highest levels of anticorrelation of ~ —0.04 predicted in the 
corresponding simply connected universe. 

As discussed earlier, the lines of sight pointing to centers of matched circles have segments of identical points and 
this can lead to anticorrelation between the CMB anisotropy in those directions. Here we explain this with the help 
of Figure ||. The segments AB' and At B of lines of sight OA and OB are identical under V. (Recall eqs.(||) and (|g) 
for the Sachs- Wolfe contribution to the correlation between temperature fluctuations in two directions in the sky.) 
The surface term depends on the correlation of the potential, $, between two points physically separated by the 
distance r±. The interference term however contains correlation between identical points, A A' and BB\ and close 
ones. The integral term has correlations between close points but is down-weighted by the extra factor of AF. Thus 
the interference term (which is usually small in simply connected models) in this particular setting can dominate the 
total CMB correlation and make it negative. Thus the existence of strongly anticorrelated spots caused by the ISW 
term is a signature of non-trivial topology. 
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FIG. 8. The negative correlation between pixels located at the center of matched circles is demonstrated. The top panel 
shows the 8 largest circle pairs for the LCH model on a sky map, as in Figure ^. The groups of pixels marked out in the lower 
sky map are distant pairs with high anticorrelation, C pp , / J C vv C v i v i < —0.3. It is visually apparent that these regions lie at 
the centers of matched circles and are labeled accordingly. 

C. Patterns Arising from Global Inhomogeneity 

The global inhomogeneity of compact spaces implies that the variance of the the gravitational potential is spatially 
dependent. If these fluctuations arose from quantum noise during inflation, $(x) would be an inhomogeneous Gaussian 
random field, with a pattern of inhomogeneity determined by the topology of the compact space. 

The CH models that we have considered also predict that the rms temperature fluctuations in the sky vary with 
direction, as is shown in Fig. || and in Fig. 10 of Paper 1. In both LCH and SCH models there are 'loud' spots - 
directions in the sky where the variance is significantly larger. For the NSW contribution, these arise when the SLS 
crosses loud regions in the Dirichlet domain where the variance £$(x, x) is large. These regions consist of points in 
the Dirichlet domain where the length of the shortest geodesic is small compared to the typical size of the shortest 
closed geodesic at other points, i.e., the closest image of the point on the universal cover is smaller than ~ R < . The 
LCH model has a loud region around a vertex of the Dirichlet domain shared by two pairs of identified faces where 
the nearest image of the points close to this vertex are at a small separation. We discussed this effect and results for 
the NSW effect in a previous paper ||. Fig. ^| shows the extent to which the ISW term masks the signature expected 
from a simple NSW consideration: the loud spot is dramatically modified, even changing sign, when ISW is included. 
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FIG. 9. The figure shows two full-sky maps of the fluctuations in the standard deviation of the predicted CMB temperature, 
plotted as pairs of 180° diameter hemispherical caps, one centered on the South Galactic Pole (SGP) and one on the North 
(NGP). The contours are linearly spaced in units of 0.03 of the mean standard deviation. The first map considers only the 
NSW contribution, where at the peak value standard deviation is 45% larger than its mean value . The second plot shows 
that the significant ISW contribution at fio = 0.6 radically alters the pattern of fluctuations in the variance since the CMB 
temperature now depends on the potential along the entire line of sight. 

There are known examples of multiply connected spaces where the pattern of inhomogeneity is a dominant observable 
effect. In [ p9| a non-compact but multiply connected horn topology was constructed and the properties of the 
perturbations were studied. It was demonstrated that there is a region in which the variance of $ is strongly 
suppressed and would correspond to a dark spot in the distribution of galaxies. For the NSW contribution to the 
CMB anisotropy, this same effect directly translates to suppression of the variance of the temperature fluctuations in 
the corresponding direction, creating a flat spot signature in the CMB anisotropy, shown explicitly for the toroidal 
horn space in ]3C| ]. 

V. BAYESIAN ANALYSIS CONSTRAINTS FROM COBE-DMR 

In this section, the goal is to explicitly evaluate how likely the various compact universe models are in light of the 
COBE data on angular anisotropies. We first review the statistical distributions for the maps derived from the COBE 
data, then show how we use our techniques to confront the data. 

The raw data of a CMB experiment comes in the form of a time stream of measurements d t at N t time-ordered 
points for each frequency channel. The data are then binned into a AT p -pixel discretization of the sky, through the 
relation dt = Pt p A p + rjt, where the Nt X N p pointing matrix P tp maps the observing time to the angular position 
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at that time and rjt is the time stream noise. A p is the true signal on the sky. With the (reasonable and checkable) 
assumption that the noise is Gaussian with covariance matrix N = (r)tT)t'), one can find the map A which maximizes 
the conditional probability V(d\A) and the pixel-pixel noise covariance matrix about it, Cat- 

A = C N P^N- 1 d, C N = (P^N^P)' 1 , (10) 

No larger moments required, given the Gaussian noise assumption. Thus, the probability of the data given the true 
sky signal A is 

where Cn — ((A — A) (A — A)^). Provided the pixelization is fine-grained enough, the A map plus the Cn contain 
all of the sky information present in the original data set dt- 

The COBE team |26| have given 6 maps at 3 frequencies, 31, 53, and 90 GHz, each of 6144 pixels of size (2.6°) 2 , 
along with information to construct C'n from the number of observations made in each pixel and the average noise 
in the radiometers over an observing time. Four our analysis, we compressed the six COBE-DMR maps p6| into a 
(A+B) (31+53+90 GHz) weighted-sum map. Galactic emission near the plane of the Galaxy sufficiently contaminates 
the primordial signal that a region ±20° from the Galactic plane is removed, along with adjacent extra pixels in which 
contaminating Galactic emission is known to be high, as advocated by the DMR team. Although one can do analysis 
with the map's (2.6°) 2 pixels, this "resolution 6" pixelization of the quadrilateralized sphere is oversampled relative 
to the COBE-DMR beam size, and there is no effective loss of information if we do further data compression by using 
"resolution 5" pixels, (5.2°) 2 p7[ ]. The celestial sphere is then represented by N p — 1536 pixels before the Galactic 
cut, with N p — 999 pixels remaining after the cut is made. Our Cn is largely diagonal, but we include the off-diagonal 
components centered on a 60° pixel-pair angle separation, which corresponds to the horn separation of the instrument. 
We remove a best-fit monopole and dipole from the cut-sky maps. Proper account is taken of the monopole and dipole 
contributions, as well as possible quadrupole contamination by Galaxy emission, by increasing the noise in associated 
template patterns pl| , p2[| . This corresponds to having arbitrary monopole, dipole and quadrupole contamination 
possible, and effectively "shorts-out" this contributions to Ct- 

Now we need a probabilistic model for the signal. There may be several physically different signals in the COBE- 
DMRdata - not only primordial CMB anisotropy, but also Galactic emission and others. However, except for the 
quadrupole contamination which we corrected for, the contribution of signals other than CMB is small in the weighted 
sum (over frequency channels) maps that we have used. We have assumed a Gaussian probability distribution for 
the primordial fluctuations of As we have seen AT/T is linearly related to $, which remains true even if all the 
effects leading to the CMB anisotropy are included. Thus, AT/T is also statistically a Gaussian random field, fully 
described by the theoretical pixel-pixel correlation matrix Ct that we have focused on in our computations: 

V{A\Ct) = {2n) NJ llCTV ^ AtC ^ < °- = C &> *») ■ M 

The Ct which enters here should have the COBE beam taken into account. We do this by calculating Ct through 
eq.(||). We used a high fc-cutoff in evaluating the C u , choosing its value to regulate the high frequency part of Ct, 
but ensured that it was much higher than the corresponding scale associated with the COBE beam size. The effect 
of the COBE beam was included by forming B^CtB, where B pp i is the COBE beam between pixels p and p 1 , related 
to the beam-shape Bi by 

B pp' = V B t p t% ■%')■ ( 13 ) 

This gives all the necessary ingredients for our analysis except for a prior probability V(Ct) for the theory, which 
may encode both our theoretical prejudices about the models and results from other observational tests. For the 
moment we leave V(Ct) unspecified and concentrate on the likelihood function of a model 

C{C T ) = V{A\C T ) = J dAV{A\A) V{A\C T ) . (14) 

The integration is carried over all virtual realizations of the sky A. The result is 

C{Ct) - (2^r^\\c N +c T \\^ • (15) 
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The likelihood, defined by eq.(14), is ultimately a function of the parameters of the model, built into the Ct- Maxi- 
mization of C(Ct) in the parameter space is a complex task, if, as is generally the case, the dependence on parameters 
is nonlinear. However, it is straightforward to determine relative likelihoods of any models with our precomputed 
C T 's. 

In the case of CH models the parameters are: the choice of manifold Ad; density parameters, of the matter O m 
and of vacuum energy, (most important is their combination Qq which sets R LS /d c , and, hence, the physical size 
of the compact space); the manifold's orientation with respect to the observed map - i.e., the triplet of Euler angles 
a; the position of the observer x. bs within the manifold; and the parameters characterizing the initial spectrum of 
fluctuations, such as the overall amplitude and the spectral tilt. If, as here, we fix the initial spectral slope, only 
the amplitude A remains free. In our studies, we have assumed a uniform prior probability for A and integrated the 
likelihood over it (i.e., marginalized the parameter). However, since C is always sharply peaked near the best-fit value 
of A, the choice of the prior for A is irrelevant. 

Thus, C(Ct) = C(Ai, fi m , Oa> ot,x bs)- The logical way to proceed would be to do many manifolds M; for each 
manifold, many different fl m and for each (M, fi m , ^a), many orientations a; etc. For the exercise presented here 
we have chosen two model spaces, one with a small volume, SCH, m004(-5,l), and one with a relatively large onw, 
LCH, v3543(2,3), and considered 'pure' open models J7a = 0. For each of these CH spaces we found the likelihood for 
three values of fl m and twenty- four different orientations. We have chosen Q m values to straddle the line i?> « R LS , 
since we have found this to be a rough boundary between models which pass and which fail the COBE-DMRtest. 

Twenty- four orientations correspond to the rotational symmetry of a cube on which the COBE— DMR pixelization 
is specified. This allows us to avoid interpolating Ct to new pixel positions for each rotation and deal only with 
remapping of the correlation matrix elements. We have not varied the position of the observer. We have chosen it 
to be at the "local maximum of injectivity radius" , from which position the space usually looks most symmetric (or 
'round'). We expect that for this observer the model will be less restricted, than for an observer at another place, 
who would see a more squashed and anisotropic space. A caveat is if the most squashed direction is partly hidden 
within the Galactic plane. 

In the torus model calculations, we chose many more manifold orientations to sample the Euler angle space, since 
Ct could be easily computed with a Fast Fourier transform using the known eigenfunctions of the torus. In general, 
we could continually refine our orientation angles to hone in on the maximum likelihood value more precisely. We 
would obviously do so if we felt that we were on the trail of a true model of the Universe, but such a refinement is 
not essential for the points we make here: that universes which are much smaller in volume than the volume within 
the last scattering surface are strongly ruled out, independent of orientation. 

In Table I we present the results for the likelihood of the compact models in our selection relative to the likelihood 
of the standard non-compact open CDM model with the same f2 m . The theory with infinite volume is known to 
fit the COBE-DMRdata well and is considered to give a good description of the data. COBE-DMRdata alone does 
not discriminate well between the infinite models with different values of f2 m and except disfavoring very low 
fi TO < 0.2, thus the choice of our reference model is not critical. 

TABLE I. The Log-likelihoods of the compact hyperbolic models relative to the infinite models with the same f2 m are 
listed. The probabilities are calculated by confronting the models with the COBE-dmr data. The values quoted are likelihoods 
marginalized over the amplitude of the initial power spectrum. The volume within the sphere of last scattering (SLS) relative to 
the volume of the compact models of the universe (V ls /Vm) is listed. The three columns of logarithm of likelihood ratios C/ Co 
correspond to the best, next best and the worst values that we have obtained amongst 24 different rotations of the compact 
space relative to the sky. The number v in brackets gives the conventional, albeit crude, translation of the probabilities to a 
Gaussian likelihood C/Cq ~ exp[— v 2 /2\. Only the last model for one specific orientation appears to be consistent with the 
COBE-DMRdata. 
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The clear conclusion to be drawn from the Table is that when Vt m is too small, so i?> < i? LS , the likelihoods are 
tiny relative to the larger models. We believe this is a robust conclusion, largely independent of the details of the 
manifold choice, orientation or the assumed spectrum of the initial fluctuations. 

What is also clear is that, near R > w R hs , it can happen that for certain manifolds and orientations, £ is higher than 
for the standard oCDM Universe. The interpretation is that some of the highly correlated spots that are predicted, 
such as those shown in Fig. ||, partly line up with the observed spots in the COBE-DMRmap. Even though there 
are realization-to-realization fluctuations, the random skies, derived from an anisotropic model with correlated spots 
built into Ct, will always be constrained to deliver pixel pairs reflecting these correlations. By contrast, in infinite 
isotropic models there are no preferred spots and in a much smaller fraction of realizations particular spot line-up 
will happen. Thus, this particular CH model at the specfic orientation would be always preferred over its isotropic 
infinite counterpart. 

Allowing the manifold and its orientation to vary, we can get this alignment from time to time, given so many 
parameters. An obvious question is how to assign the prior probability for orientation. This should obviously 
be random, except that if we actually do live in CH universe, there is a true orientation and we are not allowed 
to marginalize over orientation nor manifold choice. If high likelihood is achieved for some manifold at a specific 
orientation, one could argue that this model is a preferred explanation, at least for the COBE-dmr data. What is 
then required to test this explanation? Clearly, a strong test is to go to higher resolution. If the same manifold 
and orientation remain preferred at higher resolutions, this should spur cosmologists on to further checks of the CH 
hypothesis and search for specific signatures of the compact space. A powerful check is to search for the correlated 
circles such as in Figure ^. A manifold-independent strategy with 13 arcmin MAP data emphasized by Cornish et 



al. 1 28 1 exploits these correlated circles. A caveat is that there are other "surface terms" involving the Doppler term 
which will spoil somewhat the simplicity of this strategy. 

Figures [u] and [ll] compare theoretical realizations of the CMB anisotropy in the LCH and SCH models with 
the COBE-dmr data. They should be compared with the 'DATA' map in Fig. a Weiner-filtered picture of the 
CMB data. The Wiener-filtered map is the mean signal subject to the constraint of the observations for a theory 
characterized by a given Ct'- 

(A s \ A, C T ) = [C T (C N + Ct)- 1 ]^ . (16) 

The maps obtained with different choices for Ct in the Wiener filter often look quite similar as long as the Ct fits 
the data reasonably well. For Ct, we used that for a standard Q m — 1 CDM model, which fits the data rather well. 
As far as the visual appearance of the map is concerned, a £l m < 1 oCDM model would look very similar j3(|. The 
differences are in line with what we might expect: if there is a little more power on small scales, then the map has 
slightly more contours at small angles. When one uses a model which is greatly disfavored by the data, the Wiener 
map looks extremely different. For example, there is definitely a pronounced large-scale signal in the DMR data 
which a low fio compact hyperbolic model cannot reproduce. It then tries to interpret that large-scale signal as a 
chance (and highly unlikely) superposition of noise, which is another expression of why the model is so statistically 
disfavored. 

What one should be noting looking at the maps is the shapes of the patterns and not the specific locations of the 
patterns, since these can change from realization to realization. The full Bayesian analysis takes into account all 
possible realizations. The incompatibility of models with small Vm/V ls (SCH-^o = 0.3,0.6) is visually obvious: the 
best fit amplitudes are high which is reflected in the steeper hot and cold features. Although, the SCH-£!o = 0.9 
and LCH-fio = 0.6 models do not appear grossly inconsistent, it turns out that the intrinsic anisotropic correlation 
pattern is at odds with the data statistically (Table I). 
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FIG. 10. The figure consists of a column of three CMB sky-maps, each showing a pair of 140° diameter hemispherical 
caps, centered on the South (SGP) and North (NGP) Galactic Poles, respectively. The top map labeled DATA shows the 
COBE-dmr 53+90+31 GHz A+B data after Wiener filtering, assuming a best-fit standard CDM model, normalized to the 
COBE-dmr amplitude. The next two maps are of one random realization of the CMB anisotropy in v3543(2,3) - our choice 
of a L(arge)CH model example, for flo = 0.6 & 0.8 based on our theoretical calculations of C(q,q') convolved with the 
COBE-dmr beam. Both surface and integrated (ISW) Sachs- Wolfe effects have been included in C(q,q'). No noise was added. 
The amplitude in each model was chosen to best-match the COBE-dmr data. The theoretical sky was Wiener-filtered using the 
COBE-dmr experimental noise to facilitate visual comparison with the data. LCH with Qo = 0.8 is compatible with the data 
with a suitable choice of orientation, whereas with f2o = 0.6, it is ruled out (See Table [[]). For all the maps in Figs, [u] and 0, 
the average, dipole and quadrupole determined for this cut sky were removed. A 20° Galactic latitude cut was used, with extra 
pixel cuts to remove known regions of Galactic emission proposed by the COBE-DMRteam, accounting for the ragged edges. 
The contours are linearly spaced at 15 fiK steps. The maps have been smoothed by a 1.66° Gaussian filter. 
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FIG. 11. The three CMB sky-maps, each showing a pair of 140° diameter hemispherical caps, centered on the South (SGP) 
and North (NGP) Galactic Poles, are analogous to the lower two plots in Fig. |Io| but for the S(mall)CH model m004(-5,l). 
The fact that Qo = 0.3 & 0.6 models above are strongly ruled out by the COBE-DMRdata is obvious visually. The fio = 0.9 
model is not obviously excluded on the visual basis, but is indeed excluded on the basis of our Bayesian analysis. (See Table 

VI. CONCLUSION 

Although there are an infinite number of possible CH spaces, one can extrapolate some general conclusions on CH 
universe models from our limited exploration of COBE-dmr constraints on the LCH and SCH examples. We shall 
present the CMB constraints on a larger set of CH spaces in [[l3| . The main CMB feature of small compact universe 
models - the presence of high correlations between many well-separated pixel pairs - is also their handicap. The 
COBE-DMRdata does not, generally, favor interpretation as being a noisy random realization derived from a small 
compact model. Formally, the likelihood V (A COB e-dmr \Crpp') is much smaller for such models than for a standard 
oCDM theory in which Ct pp ' mostly just falls off with separation between the pixels. 

High correlations at large angles are numerous in a compact space with cIm < R LS and we are confident in concluding 
that such topologies are not viable models for our Universe in view of the COBE-dmr data. Of course, the possibility 
of exceptions in the infinite list of CH models remains, but the bulk of models must satisfy the, in our opinion quite 
solid, limit cLm > aR LS to pass the CMB test. Setting a — 1 would be a very conservative choice; all our numerical 
simulations are consistent with at least a = 1.4. At this limiting value the CH models we tested are excluded at 3<r 
level. 

Similar conclusions were reached by some of the authors (JRB, DP and Igor Sokolov |l^]) for flat toroidal models. 
Comparison of the full angular correlation (computed using the eigenfunction expansion) with the COBE-dmr data 
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led to a much stronger limit on the compactness of the universe than limits from other methods |33|,|34[ . The main 
result of the analysis was that Rt/R^s > 1.3 at 95% CL for the equal-sided 3-torus (with the periodicity length 2Rt, 
the diameter of the torus is cLt = 3 1 ^ 2 Rt, thus a = 2.25 ). For 3-tori with only one short dimension (or for the 
non-compact 1-torus), the constraint on the most compact dimension is not quite as strong because the features can 
be hidden in the "zone of avoidance" associated with the Galactic cut. 

Statistical properties of fluctuations in CH manifolds are anisotropic and inhomogencous. CMB predictions of 
course depend not only on the topology of the space, but also on the position of the observer and the orientation of 
the Dirichlet domain with respect to COBE-DMRsky. This lack of high symmetry is reflected, in particular, in R > and 
i?<, which are not invariant under the choice of the space basepoint (i.e., the observer position) and which quantify 
how the sky looks for a given observer. Our limit on CH models uses the invariant linear measure (Lm, but roughly 
corresponds to the condition i?> > R hs when R > is the outradius for the observer at the basepoint which maximizes 
the 'injectivity radius' of the space |l9|. This means that when the Dirichlet domain just fits into the last-scattering 
sphere, the correlation matrix is already too distorted to satisfy the data. Moving the observer to another point, will 
generally increase i?>, but will also squash the domain in some other directions. In the more anisotropic view of the 
CH space presented to such an observer, we expect to predict less favorable CMB skies. Thus, moving the observer 
away from the basepoint which maximizes the 'injectivity radius' may not relax the constraint. The absence of good 
data close to the galactic plane, (i.e., the galactic cut in the data) may help some models at specific orientations, but 
not to the extent that it does for the 1-torus space, which has an exact planar symmetry. Extensive analysis of the 
changes induced by varying observers is left to future work. 

As we emphasized in [0J^] and here, the constraints arise predominantly from predicted pattern mismatches in our 
models compared with the COBE tapestry. This is entirely encoded in Ct pp ' , which can also be expressed in terms 
of a Yg m basis. However, CTtm.i'm' is generally quite complex and reducing consideration to the isotropized Cg loses 
a substantial amount of information, since it involves a 5w projection, followed by a trace over m. More importantly, 
as we show in detail in Appendix |B], Cg has substantially increased error bars from "cosmic variance", i.e., in the 
expected theoretical fluctuations about the mean, so we can draw only extremely weak conclusions about the model. 
This is evident in the error bars on the angular power spectrum shown for a set of the CH models in Figure Thus, 
although the Cg for a compact model may fit the data reasonably well, and it sometimes does so even better than the 
corresponding infinite model with the same Oo, statistically this may be a rather meaningless observation, and, if one 
is not careful, even misleading. Some authors [ ^8|j37| have argued that because the mean Cg shape may look better 
visually, this is evidence that the models are preferred. 

To make the point quantitatively that conventional use of Cg can lead to very wrong conclusions, Table [n] compares 
the likelihood ratios for a few models obtained using just Cg information, treated as if they were statistically isotropic 
Gaussian models with the the isotropized power spectrum, with what was obtained in Table [I] when the full pattern- 
recognition statistical treatment was made. Details on the construction of the Table are given in Appendix |b| In 
all cases shown, the CH Cg is preferred over the Cg of the corresponding infinite model, but it is grossly misleading 
because of the enhanced error bars, and the huge amount of relevant information left out. All models are strongly 
ruled out, save one. And that manifold, with its specific orientation relative to the sky, is preferred even more than 
the statistically-isotropized one. 
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FIG. 12. The three panels show angular power spectra for the SCH model at fio = 0.3, 0.6 and 0.9, respectively. The larger 
error bars are the actual cosmic variance computed using the full correlation matrix information. The smaller error bars drawn 
slightly displaced to the right are the cosmic variance one would assign if one naively assumed that the Ce contained all the 
information. The excess variance reflects the incompleteness of the Ce information due to the inherent statistical anisotropy. 
The dotted (red) curve shows Ce for the simply connected 'open' universe at the same value of Qo- The error bars show the 
cosmic variance. The Ce is normalized to give the same r.m.s. power as the corresponding SCH model for the COBE-dmr beam. 
The full COBE-dmr beam, with rough Gaussian scale (T^^m ~ 17.5, has been factored out in making the plot. 



TABLE II. The Log-likelihoods of the compact hyperbolic models relative to the infinite models with the same Q m are listed. 
The probabilities are calculated by confronting the models with the COBE-DMRdata solely in terms of Ce- For easy reference 
the corresponding relative likelihood for the "best" orientation is listed from Table EL 
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Most astrophysical observations point to a matter density O m < 0.4 ||. For the matter-dominated f^o = Q m 
open cosmology, when we combine this with the dju > aR SL s COBE-dmr constraint that we sugg 
1 — tanh 2 (dx/2ad c ), we would be able to rule out topologically small CH models with d/n/d c < 2.8 (adopting the 
conservative a = 1.4). 

Recent SNIa results Jul, the emerging location of a peak in the CMB power spectrum at £ ~ 200, and the 
combination of CMB with large scale structure data jllj all point to a significant fl\ term, with less room for a small 
51o = O m + 17a, even though f2 m may be small. Ironically, although this makes open models less attractive, having a 
non-zero J7a which combines with f2 m to give Slo near unity actually returns topologically small CH spaces back to 
life, since it alleviates the constraint on the size cLm- Qualitatively, CH models with the same value of R hs /d c will 
have similar constraints, whether there is SIa or not. Some quantitative differences will of course arise because of the 
diminished contribution from the ISW effect in the £l m + Q\ model compared to the corresponding Q\ — one. 

In Figure |l3| we show the fi m — J7a parameter space, where the lines of constant R hs /d c provide rough limits on 
the viability of the CH models, depending on its size. This plot illustrates that allowing for nonzero J7a relaxes the 
limits on the allowed CH topology. This is a welcome conclusion, since one could argue that (topologically) small CH 
spaces are less complex and may be more probable for quantum processes in the early Universe to have created them, 
making them a more natural choice among other CH models p5[ . 




FIG. 13. In the Q m — parameter plane, the horizontally shaded area shows the 95% CL (confidence limit) region coming 
from the analysis of high redshift supernovae SNIa [10]. For reference, the vertically shaded region shows 95% CL restrictions 
arising in the standard infinite universes when SNIa and CMB data are combined. Since intermediate-scale CMB data dom- 
inates the CMB constraint on the first peak position and hence on f^o, where the influence of the topology would not be so 
important, the limits on Q m ,Q,A in CH spaces should be similar to those shown. We suggest that lines of constant R LS /d c 
would provide effective guides to the viability of the compact models. The two lines with enhanced weight correspond to the 
line R LS /d c = d.M/1.4 in our SCH and LCH example spaces. The allowed region is to the right of the R LS /d c = d^i/1.4 lines. 
The vertical line at fl m = 0.4 cuts out the high matter density part of the plane which is disfavored by observations. 

Although our results strongly indicate that manifolds with small R > < R LS are unlikely to survive confrontation with 
the COBE-dmr data, we emphasize that, in the i?> > R LS > i?< regime, there is some room both to have interesting 
specific CH correlation patterns and still be consistent with the COBE-dmr data. When R< is large compared to 
R LS , the results will quickly converge towards the usual infinite hyperbolic manifold results. The intermediate terrain 
still encompasses ample scope for interesting topological signatures to be discovered within the CMB. Although our 
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methods are quite general, testing all manifolds in the SnapPea census this way is rather daunting, and there are 
countably infinite manifolds not yet prescribed. What may be promising for discovery are specialized statistical 
indicators, which are less powerful discriminators than the full Bayesian approach we have used here, but not as 
manifold sensitive; e.g., the statistical techniques which exploit the high degree of correlation along circle pairs that 
p8| have emphasized, and Fig. || reveals. Maps like we have constructed will be necessary to test the statistical 
significance of such methods. We also note that dramatically increasing the resolution beyond that of COBE-dmr , to 
e.g. MAP resolution, is quite feasible with current computing power using our techniques. 

VII. ACKNOWLEDGEMENTS 

We have made extensive use of the SnapPea package and related material on compact hyperbolic spaces, which 
is available at the public website of the Geometry center at the University of Minnesota. TS acknowledges support 
during the final stages of this work from NSF grant EPS-9550487 with matching support from the state of Kansas. 
In the course of this project we have had enjoyable interactions with Janna Levin, Neil Cornish, Dave Spergel, Igor 
Sokolov, Glen Starkman and Jeff Weeks. 

APPENDIX A: INCORPORATING SMALL ANGLE CMB ANISOTROPY 

In this section, we discuss the calculation using the method of images of the primary CMB anisotropy at smaller 
angular scales where sources other than the Sachs- Wolfe effect make the dominant contribution. The most important 
effects arc the Dopplcr shift due to scattering of photons on free electrons during recombination, and a term describing 
the compression and rarefaction of photons. 

We shall consider only the scalar mode of perturbations. As in the main text of this paper, we choose to work in 
the longitudinal gauge in which the metric perturbations are described by two scalar potentials $ and W 

ds 2 = a 2 (r) [(1 + 2$) dr 2 - (1 - 2*).g.^ ) dx l dx k ] . (Al) 

This gauge is particularly suitable for the analysis of multiply connected universes, since the perturbed 3-space is 
explicitly conformal to the background 3-space, leaving the topological identification of points in the background 
coordinates x % exact. Also for the perturbed quantities in the longitudinal gauge one does not need to distinguish 
between gauge-specific and gauge-invariant definitions. 

In first-order perturbation theory, the photons observed at the origin from the direction q at the (present) moment 
To have propagated along radial null-geodesies. The position (r, x l ) on the geodesies is parametrized by (to — XtXQ) 
in terms of the affine parameter x with % — at the origin. r coincides with the present radius of the FPW horizon, 

XH- 

The transport of the temperature fluctuation Ar(r, x k , q) along the photon path is given by the Boltzmann equation 

- V X A T + (d x () A T = (d T * + q l d^) + (8 X C) (e/4 + fd^ v ) , V x = q% - d T , (A2) 

where we have used the following notations Q T> x is the total derivative along the photon path and d denotes the partial 
derivatives. Along the path the direction of the photon momentum p is opposite to the direction of q. The energy 
density fluctuations are given by local angle averaging over momentum directions e(r, x l ) = J dQpAT(T,x t ,p)/4:ir. 
The function = J* dx&Tn e is the optical depth due to Thompson scattering of photons on free electrons with a 
number density n e . The velocity is described by the velocity potential ip v . We have also omitted subdominant terms 
related to the effects of polarization and the angular anisotropy in the scattering . 

The resulting temperature fluctuation measured by the observer in the direction q is obtained by integrating eq. ( |A2| ) 
along the photon path 

At(t ,<7) = / ° D X \(d T V + q%$) e-«x) - ( e /4 + q%i> v ) d x e~^] . (A3) 



4 Other frequently used notations are v = $, tp = 5 
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Assuming the standard recombination history, the approximation of 'instant recombination' is quite accurate for 
the CMB anisotropy at large and intermediate angular scales. It assumes an instantaneous transition from the 
phase at r < t ls , where the photons were tightly coupled with the electrons, to the phase at r > t ls where they 
propagate freely after last scattering at t ls . In eq. (A2) this formally corresponds to the limit for the visibility 
function e~^ x -* — > Q(tq — t ls — x)> where 0(x) is the Heaviside step function. Correspondingly, the differential 
visibility function — d x e~^ x ^ — > S(to — r LS — x)> where 5(x) is the delta function. Omitting the monopole terms, one 
obtains the well known expression for temperature fluctuations in this limit: 



A T (r ,9) = (e/4 + $ + q l d^ v ) 



V x d T (* + *). 



(A4) 



The above expression for the CMB temperature fluctuation includes the Doppler term q l ditp v in addition to the Sachs- 
Wolfe contribution given in eq. ([|). To obtain the Sachs- Wolfe contribution in the form of eq. (|J), the relation $ = W 
valid in the (fairly good) hydrodynamic approximation to the matter content in the universe is used. Furthermore, 
the well known solution to the evolution equations for the perturbations imply e(r LS ) w — 8$(t ls )/3 in the Sachs- Wolfe 
regime. 

We have presented these well known results to stress that, in full generality, the CMB anisotropy can be expressed 
as a line of sight integral over a source function given purely in the real space. Thus, in a multiply-connected Universe 
the method of images can be directly applied to full CMB calculations in a similar manner as we have used it to 
calculate Sachs- Wolfe contribution. 

The new element that eq. (A3) has beyond the Sachs- Wolfe expression of eq. (||) is the vectorial terms in the source 
function. For the scalar perturbations all of them are spatial derivatives of the potentials along the line of sight and 
can be reduced to time derivatives through integration by parts. Indeed, q l di<& is readily converted into a 9 r $ addition 
to the integral plus the surface term <E>[t ls , (to — T LS )q], as was implicitly done both in eq. (A4) and in eq. (|^) in the 
main text. The Doppler term can also be dealt with the same way, however there is no obvious advantage of doing 
that. The reason is the (5-function like, rather than step-function like, nature of the differential visibility multiplier in 



the Doppler term. Integration by parts gives an integral term with 5r 



-C(x) 



which quickly changes sign; and having 



to integrate it numerically may not indeed be a simplification over numerical differentiation. 

As is clear in the simplifying instant recombination approximation, the Doppler term depends on the line- of- sight 
projection of a vector at last-scattering surface and not a scalar. As for a scalar field, a vector field on a multiply 
connected universe Ai u /Y, can be expressed as a sum over images of the vector field on the universal cover, M. u . 
However, in contrast to scalar fields, the discrete group T acts both on the spatial position x and the vector at that 
point, i.e., the vector at the image point 7[x] is rotated relative to the vector at x; and its projection on the line-of- 
sight is not preserved. Thus, at scales where it is dominant, the presence of the Doppler term at the last-scattering 
surface will tend to destroy the nice circle correlations that appear for the NSW case (see 
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However, the presence of vector terms is not an obstacle for the numerical implementation of the method of images. 
Inde ed, it is quit e st raightforward to implement the method of images to calculate the CMB correlation function for 
eq. (A4) or eq. (A3). For a completely general source term, the correlation function for the CMB anisotropy is the 
double integral 



C{q,q') = / d Xl 



d x 2 (S(qXi)S(q'x 2 )) 



(A5) 



In a compact (more generally, multiply connected) FRW universe, the method of images can be invoked to compute 
the source correlation function as a regularized sum over images over the source correlation function on the universal 
cover: 



{S(q Xl )S(q' X2 )r = ^MfahM-yW xM) u 



(A6) 



7er 



where the superscripts c and u refer to the quantity in the multiply connected space and its universal cover, respectively. 
Tilde refers to the possible need for rcgularization. Y is the discrete subgroup of motions which defines the multiply 
connected space and 7 fx] is the spatial point on the universal cover obtained by the action of the motion 7 G Y on 
the point x. The important point to note is that one needs to implement the action of the motion 7 on the source 
function unless all the terms in the source function are scalar quantities (as was the case for the Sachs- Wolfe effect 
that we considered in this paper) when the action is trivial. 

In summary, the CMB anisotropy correlation in a multiply connected universe can be computed in full generality 
using the method of images on the correlation function for CMB anisotropy source terms. However, it is important 
to identify the scalar, vector and tensor parts of the source functions and the action of the discrete motion has to be 
applied to non-scalar components of the source function. 
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APPENDIX B: INADEQUACY OF C e COMPARISON FOR COMPACT UNIVERSES 

In this appendix, we use the examples of the CH spaces we have studied here to discuss the limitations of a 
comparison of CMB anisotropy predictions to data solely in terms of the angular power spectrum for compact universe 
models. 

The (isotropized) angular power spectrum, Ce, widely used to summarize the CMB anisotropy predictions of a 
theoretical model, is defined by 



£{£ + !) 



2tt(2£ - 



) - . , f AT 

Yy ^2 {a tm a* lm ), where a im = j dQ, s —(q)Y im (q) . (Bl) 



The ( ) denotes an ensemble average of the random variable enclosed. In full generality, the expectation values of the 
pair products of the spherical harmonic coefficients (ae m a^, m ,) are related to the correlation function C(q, q') by 

(a em a* e , m ,) = JdSl 9 J dSly (^f(q)^f(q')^ Y im (q)Yj m ,tf) = J dSl 9 J dfy C{q,q')Y lm {q)Yl, m ,{q') . (B2) 

In comparing theoretical predictions to observations, the power spectrum, Ce, is a useful "data compression" of 
the N(N + l)/2 numbers in the N x N correlation matrix, Ct, into » N multipoles values of Ce . However, it 
is sometimes overlooked (see e.g., fl38| , |37| ] ) that a necessary condition for this compression to be "loss free" is that 
the CMB fluctuations are statistically isotropic, i.e., the ensemble averages such as C(q, q') or (ae m &em) are invariant 
under rotation. 

The standard (simply connected, FRW) universe models respect global isotropy and predict a statistically isotropic 
CMB sky. The full angular correlation function C(q, q') between two directions is then solely a function of their 
separation, q ■ q' . Using eq.(B2), the equivalent statement is that (ae m a}, m ,) = 2nCe / (£{£ + 1)) See' 5 mm ', is diagonal in 
the (£, m) space and independent of m. The angular power spectrum Ce is uniquely related to the correlation function 
through the expansion (^) and contains the same information. 

All Euclidean or Hyperbolic compact spaces violate global isotropy and the CMB temperature fluctuations are 
statistically anisotropic, i.e., C'(q, q') ^ C(q ■ q'). As a consequence, the (aemO,}, m ,) now define the relevant angular 
pow er s pectra, and these are are no longer required to be diagonal and independent of m. Thus, Ce as defined by 
cq. ( |Bl[ ) misses both off-diagonal pair products of a^ m 's and isotropizes further by summing up the m-dependent 
diagonal ae m pair products. 

Figure shows the normalized full angular power spectrum 



'm')/ V {aim a* tm ){av m < a}, m ,) (B3) 



Pirn = \ a im a p 

for all the (£, m) with 2 < I < 10 for the CH models considered in our paper. This cross correlation coefficient 
clearly shows the the existence of significant off-diagonal (aem o|' m ') m the m004(-5,l) at £Iq — 0.3 where p\™ ranges 
between —0.5 up to 0.6. This is the most extreme case in our set of models since it has the smallest size of the SLS 
relative to the Dirichlct domain (Vls/^m = 153.4) amongst the models discussed in this paper. As Vls/Vm decreases, 
the number of (I, m) pairs with high p\™ decreases and the (ae m a *e'm') matrix becomes closer to diagonal. However, 
for the same CH space the range of p\™ remains comparable, e.g., for m004(-5,l) at Oq = 0.9 (Vls/^m = 1-2), Pem 
ranges from —0.42 to 0.5 for 2 < i < 10. The larger CH v3543(2,3) model is more isotropic, the (ae m ag, m ,) are 
somewhat more diagonal with most of the significant off-diagonal terms at low I. The ranges in p\™ are ±0.29 and 
±0.23 for flo = 0.6 and 0.8, respectively. The m-dependence of the diagonal products (aema*e m ) is shown in Figure 111. 
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SCH: m004(-5,l> [£2 = 0.300] 
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15 25 35 45 55 65 75 85 95 105 115 




95 105 115 



l(l+l)+m+l 



l(l+l)+m+l 



FIG. 14. The figure demonstrates the nondiagonal nature of the expectation values of ai m pair products when the CMB 
anisotropy is statistically anisotropic using the examples of the SCH model at £Iq = 0.3 (left panel) and £1q = 0.9 (right 
panel). We use a single positive integer index, n = £(£ + 1) + m + 1 to uniquely represent each (l,m). Each square 
represents a pair of (£, m) spherical harmonic indices for all (£, m) for 2 < £ < 10. The cross correlation coefficient, 
p\™ = (a lm d\, m ,) I \l (ai m a^ m )(aii m i a*, m ,) at a given (n, n') square is represented by the level of grey (or color) shown 

in the accompanying palette. (The squares corresponding to pairs with mild cross correlation, \p\^ \ < 0.12, have been left 
blank to highlight the strong ones) . For statistically isotropic CMB anisotropy only the diagonal p^™ = 1 terms will be nonzero. 
As expected, the CMB anisotropy in the SCH model is more isotropic at Hq = 0.9 as evident from the more diagonal nature 
of pf™ in the right panel. 
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FIG. 15. The three panels of the figure show the relative deviation of 1(1 + l){aen 



,) from the mean value Ct for all 



(£, m) ( 2 < £ < 20) in the SCH model at f2o = 0.3, 0.6 & 0.9, respectively. The non zero deviations demonstrate that the 
diagonal (at m a^, m ,) are not independent of m when the CMB is statistically anisotropic. The relative deviations decrease with 
fi m as the space become larger relative to the SLS. The point with large deviation in the upper two panels corresponds to 
(£,m) = (8,0). 
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For statistically anisotropic CMB fluctuations, the Cg contains less information than the full correlation matrix 
C(q, q') independent of the underlying statistics. For concreteness, we now focus our discussion on Gaussian random 
CMB fluctuations which are completely specified by C(q,q'). For statistically anisotropic CMB, the incompleteness 
of the information contained in the Ct is reflected in the enhanced cosmic variance, (Cf ). Let us s plh the correlation 
matrix into an isotropic part determined by Cg through eq.(Q) and an anisotropic term containing the remainder, i.e., 



C(q,q') = C 1 (q i q') + C A (q,q') 
so that, by definition, the isotropic part C l (q, q') is described by the set of coefficients Cj of Legendre series 



C\qA') = Y. 



£ + 1/2 
t(l+l) 



C e Pi(q ■ q') , 



(B4) 



(B5) 



and the anisotropic part C A (q, q') is orthogonal to the Legendre polynomials 



<m 9 / d^C A {q,q')P t {q-q') = 0. 



(B6) 

The presence of the nonzero anisotropic part C A (q, q') is the main attribute of statistically anisotropic models, and 
in particular of the CH models. Consider now the effect this term has on the probability distribution P(Cg) of the 
isotropized power spectrum estimator 



C 



5 _ + 



) f f AT AT 

I) E j J d^. — {q) — {q')Y lm {q)Yl m {q'). 
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Of course, the expectation value of the estimator is determined solely by C l (q, q') (as it should be using eq.(p(i|)) 



(B7) 



<6> = ^±^ / <n\„ j ,m, l .c{ q .< 1 ')P l M-<h- 



(B8) 



However, the distribution P(Ce) depends on the anisotropic part C A (q,q'). In particular, the variance of Ci, which 
can be calculated as a four-point correlation of the AT's, is enhanced due to the influence of C A (q, q'): 



r(C<) = (Cj) - (Ct) 2 = 2 
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(B9) 



In the final expression, the first term is the well known result for the cosmic variance of Ce, strictly valid only for 
statistically isotropic CMB fluctuations. The second term was obtained using the fact that C A (q\, q 2 ) and Pe(qi ■ q 2 ) 
are symmetric functions. It represents a positive definite correction to the standard cosmic variance arising from the 
anisotropy of C(q,q'). Hence, the cosmic variance is always larger for a statistically anisotropic CMB compared to 
that expected for the statistically isotropic case with the same Cg_. 

To calculate cq.(B9) numerically, we took our Ct pp > and used matrix representations of the Pi at the relevant 
COBE-pixel pairs. As a check of accuracy, the same procedure was applied to the infinite statistically isotropic 
models and the first term in eq.(B9) was recovered. This is a stringent test of the cancellations required to obtain 
this result. We found that for the finite compact models as well as the infinite models, our results are accurate for 
angular scales bigger than the beam, including the finite pixel effects. 

Figure |l2| shows the angular power spectrum for a set of the CH models with the associated cosmic variance. It 
demonstrates that the cosmic variance of Cg in CH models is significantly higher than one may naively assign assuming 
statistical isotropy. As a result, Ci& do not strongly distinguish CH models from the corresponding standard isotropic 



5 For simplicity we ignore the experimental window function, W(q, q'). For COBE— dmr it is isotropic, with We = Bf, (where 
Bi is spherical transform of the isotropic experimental beam,) and the effect of including the window function is to scale the 
Ct by Bj. 
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models and, on their own, are not very restrictive. Comparing CMB predictions of CH spaces to data using Ce alone 
is not incorrect, but inadequate since one has not used the full information available. The larger cosmic variance 
simply implies that the theoretical prediction is weaker. The argument that a comparison based on more information 
is more relevant is quite obvious. Any evaluation of the relative likelihood of a model based on Cg is superseded by 
a likelihood analysis that uses the complete correlation information. A model that does well solely in terms of Ci but 
fails in terms of the correlation matrix Ct pp ' should be considered ruled out. 

In Table || we compare the relative likelihood for our models obtained using Ce versus that obtained using the full 
correlation information. 

The procedure we use to determine the first line in the Table is to take the compact model's Ct pp > , use the P( 
matrices to calculate the theoretical Cg, as in eq.( p38| ), then assume that we have an infinite model with exactly that 
power spectrum, so that these CgS encode all the information in the theory. We then use these to calculate new Ct pp < 
for this theory, and determine the full Bayesian likelihood for it relative to the true infinite open model with the 
same fio and its corresponding Ci. The results show that at roughly the +la level, these power spectra are preferred. 
Note that the likelihood ratios in row 1 are independent of orientation. However, the second row of the Table shows 
that when the full information on angular patterns is included in the analysis, the likelihoods change dramatically, 
strongly disfavouring small compact models, as in Table [| As well, the case for the Qq = 0.8 LCH model at the best 
orientation is enhanced by the inclusion of the full pattern information. 

Models which fare very poorly with respect to a full correlation comparison may well look favored based on the 
C{. The reason for this is not hard to understand. If the compact space is not much larger than the SLS then it 



predicts strong anisotropic correlation features in the CMB sky which are at odds with the data (see § IV). However, 
an isotropized measure such as the Ce is insensitive to these features. This implies that the comparison of CMB 
anisotropy in CH models using Cg alone is grossly inadequate and could be quite misleading. 
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